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A detailed investigation of the unusual dynamics of the magnetization of (Et4N)3Fe2Fg (Fe2), 
containing isolated [Fe2Fg] 3_ dimers, is presented and discussed. Fe2 possesses an S = 5 ground 
state with an energy barrier of 2.40 K due to an axial anisotropy. Poor thermal contact between 
sample and bath leads to a phonon bottleneck situation, giving rise to butterfly-shaped hysteresis 
loops below 5 K concomitant with slow decay of the magnetization for magnetic fields H z applied 
along the Fe-Fe axis. The butterfly curves are reproduced using a microscopic model based on the 
interaction of the spins with resonant phonons. The phonon bottleneck allows for the observation of 
resonant quantum tunneling of the magnetization at 1.8 K, far above the blocking temperature for 
spin-phonon relaxation. The latter relaxation is probed by AC magnetic susceptibility experiments 
at various temperatures and bias fields Hoc- At Hoc = 0, no out-of-phase signal is detected, 
indicating that at T > 1.8 K Fe2 does not behave as a single-molecule magnet. At Hoc = 1 kG, 
relaxation is observed, occurring over the barrier of the thermally accessible S = 4 first excited state 
that forms a combined system with the 5 = 5 state. 



PACS numbers: 75.45.+j, 75.50.Xx, 75.60.Ej, 75.60.Jk 

I. INTRODUCTION 

Among the most fascinating aspects of magnetism in 
recent years has been the discovery that the magnetic 
moment of individual molecules can give rise to magnetic 
hysteresis phenomena. Such systems lie at the interface 
between classical and quantum-mechanical regimes, and 
are thus of great fundamental as well as practical in- 
terest regarding future applications related to quantum 
computing^ 

Molecular clusters such as Mni 2 ) 2 ' 3 ' 4 ' 5 Fe$£ and Mn^£ 
have been found to exhibit hysteresis curves character- 
ized by a remanence. They have been dubbed single- 
molecule magnets (SMMs) as each molecule behaves as 
a single-domain nanomagnetic particle. These high-spin 
systems possess an energy barrier A = l-DsIS 2 for re- 
versal of the direction of the magnetic moment ( "spin" ) , 
arising from the combination of a large ground state to- 
tal spin quantum number S with a uniaxial anisotropy 
(negative axial zero-field splitting parameter Ds). At 
kT <C A the magnetization can be blocked and relaxes 
very slowly, giving rise to the hysteresis. Many SMMs 
such as Mni^*4 and Feg£ exhibit quantum tunneling of 
the magnetization (QTM) through the energy barrier. 

Alternatively, a different type of hysteresis lacking a re- 
manence has been observed for [Fe(salen)Cl]2 8 and more 
recently for the S =1/2 spin cluster VigS as well as the 
antiferromagnetically coupled (S = 0) ferric wheels FegiS 
and Fei2>ii In these low-spin systems the ground state 
does not possess a magneto-structural energy barrier for 
the reversal of the total spin>i£ Instead, the hysteresis 
is created by dynamically driving the spin system out 
of thermal equilibrium due to insufficient thermal con- 
tact of the sample with the heat bath, thereby causing a 
phonon bottleneck (PB) effect. It possesses a characteris- 
tic butterfly shape that arises from a fast spin reversal at 



the field value where the anticrossing of the two lowest- 
energy levels occursi2ii& Generally, the magnitude of the 
anticrossing level splitting determines the quantum dy- 
namics of the spins. In the low-spin systems, this split- 
ting is relatively large (>10~ 3 K), allowing for fast spin 
reversal yielding the butterfly shape. On the other hand, 
in high-spin systems the tunnel splitting between the low- 
est levels of typically < 10~ 6 K renders QTM too slow for 
allowing a butterfly-shaped hysteresis. Therefore it was 
previously argued that butterfly hystereses are not ex- 
pected for high-spin systems^ Nevertheless, recently for 
the SMM Mng a butterfly hysteresis caused by the PB 
effect was reported above the blocking temperature for 
SMM behavior, but no analysis was offered^ Very little 
is known to date about the physics of this phenomenon 
in spin systems with an energy barrier. However, such 
information is of particular interest as the presence of 
the barrier principally enables these clusters to behave 
as single-molecule magnets. 

The tri-/i-fluoro bridged [Fe2Fg] 3 ~ dimer molecule in 
the compound (Et4N)3Fe2Fg (Et4N = tetraethylammo- 
nium), abbreviated Fe2, possesses an S = 5 ground 
state with an energy barrier of 2.40 K for spin rever- 
sal (Fig.QJ As reported previously^ Fe2 exhibits slow 
relaxation of the magnetization and magnetic hysteresis 
below 5 K. The hysteresis adopts a butterfly shape tied 
up at H z = despite the presence of the energy barrier. 
The observation of QTM prompted us to ascribe these 
phenomena to the unusually slow relaxation of a single- 
molecule magnet^ However, meanwhile additional data 
have provided compelling evidence that the observed re- 
laxation depends on the thermal insulation of the sample. 

Here we present a detailed study on the magnetic 
properties of Fe2. The earlier findings are reinterpreted 
in terms of the phonon bottleneck. A microscopic 
model based on the rapid absorption/emission of reso- 
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FIG. 1: Double-well potential of the 5 = 5 ground state of 
Fe 2 with the energy barrier A = \D S S 2 \ = 2.40 K. Full and 
broken arrows indicate spin-phonon and QTM transitions, re- 
spectively. 

riant phonons is developed that allows for an accurate 
reproduction of the observed hysteresis curves. Finally, 
alternating-current (AC) magnetic susceptibility data re- 
flecting the fast relaxation behavior of Fe2 unaffected by 
the phonon bottleneck effect are presented and discussed. 



II. EXPERIMENTAL 

Needle-shaped single crystals of (Et4N) 3 Fe 2 F9 (Fe 2 ) 
up to 5 mm in length were prepared according to Ref . Il4l 
The dimer symmetry is exactly C^h with the threefold 
axis lying parallel to the hexagonal crystal axis cMl 
The intradimer Fe'Fe distance is 2.907 A. The dimer 
molecules are well separated from each other, leading to 
interdimer Fe'Fe distances of at least 8 A, thus making 
interdimer interactions extremely inefficient. 

Samples for magnetic experiments were prepared as 
follows: To prevent hydrolyzation, samples were pre- 
pared in a drybox under a nitrogen atmosphere. All the 
crystals were perfectly transparent, indicating that no 
hydrolysis occurred. A single crystal of 2.46 mg (sample 
A) was sealed in a glass tube under N 2 in H || c (H z ) ori- 
entation. At liquid helium temperatures, the N 2 is frozen 
out and its vapor pressure inside the glass tube becomes 
very low, leading to a poor heat contact of the crystal 
with its surroundings. Alternatively, three single crys- 
tals (total mass 48 mg) were embedded in Apiezon grease 
in a dried gelatine capsule to yield a sample (B) with 
improved heat contact. Alternating-current (AC) mag- 
netic susceptibility measurements were performed both 
on sample B and on a polycrystalline powder sample (C) 
of 25 mg made of crushed single crystals, sealed in a dried 
gelatine capsule. 

Magnetic measurements were carried out using 
a Quantum Design MPMS-XL-5 Superconducting 
Quantum Interference Device (SQUID) magnetometer 
equipped with a 5 T magnet. The oscillation frequency 
of the AC field was varied between 10 and 1488 Hz at an 
AC field strength of 1 G. All data were corrected for the 
diamagnetism of the atoms using Pascal's constants. 



FIG. 2: Dependence of the relaxation time T2 on H z at T = 
1.8 K with error (3a). The dotted line is a guide to the eyes. 



III. RESULTS 

A. Relaxation of the Magnetization 

Magnetization relaxation experiments were performed 
using sample A in H || c orientation at fixed values of 
H z and bath temperature T, after saturation at 30 kG 
and subsequent quick (« 150 G/s) reduction to the in- 
dicated H z value. Relaxation times t 2 were extracted 
from least-squares fits of a single-exponential law to the 
relaxation curves for various values of H z . The resulting 
field dependence of r 2 at 1.8 K (Fig. EJ reveals that r 2 
varies between 240 s at H z — and 870 s at 3.3 kG. 
Importantly, 3 distinctive dips are observed at H z = 0, 
1.1, and 2.7 kG. 



B. Butterfly Hysteresis Curves 

Magnetization measurements at fixed T were per- 
formed in H || c (H z ) orientation by varying H z from 
+50 kG to -50 kG and back to +50 kG. For the crystal 
embedded in Apiezon grease (sample B), the magnetiza- 
tion follows the thermal equilibrium curve (asterisks in 
Fig. |3J). Alternatively, for a crystal mounted in the glass 
tube (sample A), it features a butterfly-shaped hystere- 
sis loop characterized by two distinct features: (i) the 
magnetization at H z = is zero, and (ii) near H z = 
the magnitudes of both the magnetization detected while 
approaching H z = and the one while receding from 
it are larger than the thermal equilibrium value (inset 
of Fig. |3J. As this butterfly hysteresis is perfectly cen- 
trosymmetric, further data in Fig. 0] show only the posi- 
tive wings. The hysteresis effect decreases with increas- 
ing T, vanishing above 5 K (Fig. 0J. Reduction of the 
sweeping rate from r k 3 G/s to «1 G/ai£ narrows the 
breadth of the wings (Fig. 2 in Ref. fl5f) . suggesting that 
for a substantially smaller value of T the magnetization 
would follow the thermal equilibrium curve (dashed line 
in Fig. 0). 

A drastic change in the magnetization curve of sample 
A occurs upon a very fast sweep of the magnetic field. 
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FIG. 3: Hysteresis of the magnetization of Fe2 single crys- 
tals (open circles: sample A, asterisks: sample B) at 1.8 
K obtained with V ~ 3 G/s. The arrows indicate the di- 
rection of the measurement. The solid lines are calculated 
using the resonant-phonon model (Sec. II V CI . The dashed 
line represents the thermal equilibrium curve calculated with 
J = -2.23 K, D = -0.215 K, and g = 2.00. 14 The inset 
highlights the region close to H z — 0. 
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FIG. 5: Magnetization data (sample A) at 1.8 K obtained as 
follows: after saturation at +30 kG the field was quickly (r ~ 
150 G/s) reduced to the given H z values including negative 
ones, and the magnetization was measured within 50 s after 
reaching H z . The solid line is the adiabatic curve calculated 
using the resonant-phonon model described in the text with 
H c ff :X = 350 G, while the dashed line represents the thermal 
equilibrium curve. 
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FIG. 4: Positive wings of the butterfly hysteresis of the mag- 
netization of a single crystal (sample A) for H || c orientation 
at different temperatures as indicated with F w 3 G/s. The 
solid lines are calculated using the resonant-phonon model 
(Sec. WU . 



Fig. El shows magnetization data obtained as follows: af- 
ter saturation at +30 kG, the field was quickly (r « 150 
G/s) reduced to the given values of H z including nega- 
tive ones, and the magnetization was measured within 50 
s after reaching H z . During such a fast sweep, the mag- 
netization essentially switches from positive to negative 
saturation upon crossing H z = 0. 



C. AC Magnetic Susceptibility 

To complement our studies on the magnetization re- 
laxation in Fe2, alternating-current (AC) susceptibility 
experiments were performed on both the polycrystalline 
sample C and the single crystals in Hac II Hdc II c orien- 
tation embedded in Apiezon grease (sample B). -Hdc de- 
notes a static bias field and for sample B corresponds to 
H z . With respect to relaxation, no significant differences 
were found between the two samples. 18 Experiments were 
conducted within a temperature range 1.8 < T < 15 K 
with an amplitude of -Hac = 1 G and oscillation frequen- 
cies 10 < v < 1488 Hz. Importantly, an out-of-phase 
susceptibility (x") was observed only upon application 
of an additional static magnetic field -Hdc ■ 

Figure [S] shows the frequency dependence of x" v s T 
at selected static fields of 1 kG and 5 kG using the poly- 
crystalline sample. As v is reduced from 1488 Hz to 100 
Hz, at 1 kG the peak maximum T max shifts from 5.2 K 
to 2.45 K, and further reduction to 10 Hz causes it to 
shift below 1.8 K. At 5 kG the situation is similar. In 
both cases Tmax increases with v, from which the tem- 
perature dependence of the relaxation rate 1/r = iiru 
is derived, shown in Fig. 0as Arrhenius plots of ln(l/r) 
versus 1/T max . At Hdc = 1 kG the data follow the Ar- 
rhenius law 

where l/r = 1.0 x 10 5 s _1 is the intrinsic relaxation rate 
and AE — 12.6 ± 0.2 K the kinetic energy barrier. At 
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FIG. 6: Out-of-phase AC susceptibility x" vs T at Hoc = 
1 kG (top) and 5 kG (bottom) for a polycrystalline sample 
(C) of Fe2 in a 1 G AC field oscillating with the indicated 
frequencies v. 



Hue — 5 kG, however, the temperature dependence of 
1/t is linear (inset of Fig. 0). 

Interestingly, for higher frequencies, x" is considerably 
larger at 5 kG than at 1 kG. Indeed, we observed that 
for v — 997 Hz, x" increases linearly with Hue below 2.5 
kG and appears to level off at higher fields (not shown). 



IV. ANALYSIS AND DISCUSSION 
A. Phonon Bottleneck Effect 

In H z orientation, Fe2 single crystals in a glass tube 
(sample A) exhibit slow magnetization relaxation on the 
timescale of 10 2 — 10 3 s, whereas crystals embedded in 
Apiezon grease (sample B) do not. Obviously this relax- 
ation depends on the degree of thermal insulation of the 
sample. Poor thermal contact between the crystal and 
the heat bath (as for sample A) is well known to lead 
to the macroscopic observation of the phonon bottleneck 
(PB)i 19 i 20 Also, the butterfly hysteresis loops observed 
for Fe2 closely resemble those recently reported for the 
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FIG. 7: Temperature dependence of the AC x" relaxation rate 
1/t = 2-kv shown as an Arrhenius plot (In 1/t vs 1/T max ) as 
extracted from the data in Fig. |S] The solid line represents 
a least-squares fit of the Hoc = 1 kG data to the Arrhenius 
law (Eq.0, yielding the prefactor 1/to and the kinetic energy 
barrier AE as indicated. The inset shows the same data on 
a linear scale. Note the linear dependence of 1/t on T max for 
the Hoc = 5 kG data. 

spin clusters Vi5j£ Fe6r^ an d Fei^ii that have been at- 
tributed to this effect. Thus, it is the phonon bottleneck 
that causes the observed slow relaxation in Fe 2 . 

The energy exchange between the spin system and the 
bath occurs via the phonons in the crystal. The PB 
describes the fact that at low temperatures the num- 
ber of spins is much larger than the number of available 
phonons. Thcrmodynamically speaking, the heat capac- 
ity of the spins C s far exceeds that of the phonons C v *^ 

b=^±^^»X , (2) 

such that b w 10 4 — 10 6 i2i The observable relaxation time 
is2i 

T2 = T sp + &T pb , (3) 

with T sp and T p b denoting the spin-phonon and phonon- 
bath relaxation times, respectively. Importantly, upon 
poor thermal contact between the phonons and the bath, 
Tpb becomes very large such that T2 ~ br p h , thus render- 
ing the bottleneck macroscopically observable. In a PB 
situation, the initial rapid energy transfer from the spins 
to the phonons quickly heats the latter to the tempera- 
ture of the former, such that these phonons become in 
resonance with the energy differences between individual 
\M) levelsi^i Thus, the spins and these resonant phonons 
form a single coupled system that can only very slowly 
exchange energy with the bath, leading to the observed 
magnetization relaxation. 
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B. Quantum Tunneling of the Magnetization in Fe2 

A very important result of our studies is the observa- 
tion of dips in the field dependence of r 2 (Fig. |2J) . Such 
dips are often observed in single-molecule magnets such 
as Mni2j2^ and interpreted as fingerprints for resonant 
QTM at the respective field values, which lowers the re- 
laxation time. In the case of Fe2, the dependence of T2 on 
H z appears very intriguing as the slow relaxation arises 
from the weak thermal contact between the sample and 
the bath, which hardly depends on H z . However, since a 
tunneling pathway represents a bypass, QTM occurs in 
parallel to thermal relaxation involving phonons. 23 Thus, 
its presence can significantly alter the observed relaxation 
time T2 even if the latter is dominated by the phonon- 
bath relaxation. Indeed, recently Chiorescu et al. re- 
ported magnetization relaxation data on the V15 spin 
cluster under virtually adiabatic conditions analogous to 
our situation for Fe2*& The authors found that around 
H z = the measured phonon-bath relaxation time T2 
dropped by a factor of 2 — 3 due to the level anticrossing. 

In the C^h dimer symmetry of Fe2r^the full anisotropy 
Hamiltonian including higher-order terms reads^i 



s 2 



1) 



(4) 



305(5 + 1)5 2 2 + 255 2 - 65(5 + 1) 
231$; - 3155(5 + 1)5 2 4 + 7355* 
5255(5 + 1)5 2 + 2945 2 - 55 3 (5 



where Ol = 355* 
35 2 (5 + l) 2 , 61 
1055 2 (5 + 1) 2 5 2 
l) 3 + 405 2 (5+l) 2 - 605(5+1), and 6% = (5+ + S ( L)/2. 
The BlOl term mixes wavefunctions with AM = ±6, 
allowing for resonant tunneling at H z = between the 
\M) levels -3/+3 (Fig. [IJ and at applied fields when the 
-2/+4 and -1/+5 levels cross. For B® = Bq = 0, the 
fields H Z IM at which resonant tunneling is expected for 
the 5 = 5 state in Fe2 are given byS 3 - 



H 



MM' 



nD 5 



= n ■ 714 G 



(5) 



where n = M+M' is an integer ranging from to 9. Res- 
onant tunneling may therefore occur for n = 0, 2, and 4 
at field values of H z = 0, 1.4, and 2.85 kG, correspond- 
ing to the \M) level crossings -3/+3 (Fig. [J), -2/+4, and 
-1/+5, respectively. Therefore we can expect three dips 
in the field dependence of T2, in agreement with the ex- 
perimental data in Fig. |3 The dips are centered at 0, 
ral.l, and «2.7 kG instead of 0, 1.4, and 2.85 kG, sug- 
gesting nonzero values for and Bq in Eq. (0J. Thus 
the observed dips in Fig. [2] are consistent with quantum 
tunneling processes in Fe2 at 1.8 K. However, QTM does 
not appear to be dominant in Fe2 as the drop of ti due 
to QTM compared to pure thermal relaxation at nonres- 
onant fields is only on the order of 25 % for H z = 1.1 and 
2.7 kG and 75 % for H z = (Fig. EJ) . This behavior is in 
contrast to the situation in, e.g., Mni2 in absence of a PB 
effect, where upon QTM the relaxation time decreases by 



several orders of magnitude i 3 ^ We ascribe this difference 
to the fact that in Fe2 the tunneling is phonon-assisted 
and thus slowed down by the bottleneck effect. 

Importantly, our results on Fe2 demonstrate that quan- 
tum tunneling of the magnetization can be observed even 
though the relaxation monitored mainly represents the 
phonon-bath relaxation. By creating a PB situation, 
QTM is observed at 1.8 K, i.e. far above the blocking 
temperature for magnetization relaxation reflecting in- 
trinsic spin-phonon relaxation. 



C. Modeling the Butterfly Hysteresis Loops 

The butterfly hysteresis reflects a spin reversal at 
H z = 0. For [Fe(salen)Cl]2, it was reproduced using a 
phenomenological, thermodynamic modelj^ whereas for 
Vig^ and FegiS a dissipative two-level model with a level 
anticrossing was employed. In contrast to these systems, 
for Fe2 the presence of 11 levels and the energy barrier 
for spin reversal in the 5 = 5 state complicate matters 
considerably. As QTM is significant but not dominant 
in Fe2 (Sec. II V B|) . it is likely that the spin reversal is 
predominantly achieved by thermal activation over the 
barrier. To reproduce the butterfly shape of the hystere- 
sis curves in Fe2, we thus decided to employ a novel, mi- 
croscopic model that explicitly accounts for the barrier. 
Note that although this model is based on the absorption 
and emission of resonant phonons, it does not exclude the 
possibility of QTM. 

As the resonant phonons cannot relax with the phonon 
bath to regain thermal equilibrium, they produce coher- 
ent transitions of the spins on both sides of the barrier 
at all values of H z ; however, only at H z = these tran- 
sitions lead to a change in the macroscopically observed 
magnetization. In the adiabatic limit the total energy 
of the combined system remains constant. Therefore the 
number of the absorptions and emissions has to be equal, 
leading to the spin flip at H z = 0. The coherent spin 
transitions AM = ±1 occur between all the spin lev- 
els \M), M = — 5, . . . , 5. They can be mimicked math- 
ematically by effective oscillating transversal magnetic 
fields H e s,x[t) = H cS . x J2t 5 5 cos(u> M ,M+it). In its sim- 
plest form, the resulting effective Hamiltonian for the 
5 = 5 state of Fe2 in H II c orientation reads 2 ^ 



H c ff = DS Z + gnBH e s,x(t)S x + g^ B H z S z 



(6) 



Following the derivation of the Hamiltonian in the gener- 
alized rotating framed (details see Appendix) , we obtain 



ay grot 
"■off 



gf-iBH c ff. x S x + g^BH z S z 



(7) 



The levels \m) = —5, ... ,5 correspond to the combined 
spin+phonons system. Their response ?i^ ot is shown in 
Fig. El for H eS , x = 350 G. Within about -0.5 < H z < 
+0.5 kG, they become strongly mixed, and at H z = 
they are equidistantly split by lOg/is-f/cff.x = 0.47 K. 
Note that in the generalized rotating frame the energy 
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FIG. 8: Dependence of the 11 energy levels \m) denned 
by 7-^g 0t (Eq. on H z , arising from the resonant phonons 
(^ e ff,x = 350 G). The inset highlights the anticrossing at 
H z =0. 



of a given resonant phonon is added to the energy of 
a given \M) levels Consequently, the energy difference 
to the adjacent level and hence the entire energy bar- 
rier disappear. Therefore, the resonant phonons provide 
an efficient way to overcome this barrier, leading to the 
butterfly-shaped hysteresis. 

In the adiabatic limit corresponding to the data in 
Fig. 03 obtained upon a fast (r » 150 G/s) sweep, the 
combined spin+phonons system remains decoupled from 
the bath. Thus the magnetization essentially switches 
from positive to negative saturation upon crossing H z — 
0, indicating that during the entire sweep only the low- 
est level in Fig. [H] is populated. Consequently, M(H Z ) oc 
ds mmsa /dH z , where £ mmax is the energy of the lowest- 
energy level defined by Wgf (Fig.©. With iT efi> = 350 
G the agreement with the experimental data is excel- 
lent (Fig. |SJ, thus lending credence to our model, in 
which i?off ,x is the only adjustable parameter. The range 
—0.5 < H z < +0.5 kG in which the spin reversal occurs 
(Fig. [SJ mirrors the zone in which the levels are highly 
mixed (Fig. [SJ . Thus it directly depends on the magni- 
tude of H cS , x by {dM/dH z ) H ^ oc l/H eSiX . 

In nonadiabatic situations with a slow sweeping rate, 
corresponding to the butterfly hystereses in Figs. |3] and 

the combined system is able to partially relax to ther- 
mal equilibrium during the measurement. Thus, the fast 
reversal process has to be combined with the slow relax- 
ation process. Inside a region around H z = we apply 
the resonant phonon model (vide supra) with energy lev- 
els e m of the combined system to account for the fast 
spin reversal. Additionally, the slow sweeping rate to- 
gether with the weak phonon-bath coupling gives rise to 
the presence of a small number of nonresonant or thermal 
phonons, which allows for the slow relaxation of the com- 
bined system to the bath. Thus, at H z — the relaxation 



time becomes 23 



1 



2S ( £ 

E- 



m j+1 



-, (8) 



where m m i n = mi and m max = m2s+\- The spin-thermal 
phonon transition rates read 



W mj+1 



9oSmj + (£ 



m j+1 



48npv 5 h 4 e (^ j+1 -e mj )/kT - 1 

2qm j+1 / \3 (e mj ~e m )/2kT 

9o D mj V fc "ij + l t m j ) e 



24Trpv 5 h 4 sinh[(e 



V2kT] 



(9) 



where go is the spin-thermal phonon interaction parame- 
ter, p = 1.36 x 10 3 g/cm 3 the mass density^ v the sound 



velocity, and S„ 



(S - Mj)(S + Mj + l)(2Mj + l) 2 



in the quantization axis defined by H c s yX (note that 
H z = 0)m& Considering only transitions between the two 
lowest levels, Tj n can be approximated by 



1 



sinh[(e„ lmax - e m2S )/2kT] e^s" 



c )/2feT 



1 + 



J/kT 



(10) 



5o^m?s aX ( £ r 



e m2S ) 3 /247rpw 5 ^ 4 = 



where 7 ™»- 

810gg (gn B H e s, x y /2'inpv 5 h' i at H z = 0M. Notably, as 
for a two-level system m,2s — m-mim this equation shows 
the well-known tanh behavior for the relaxation time of 
a two-level systemiSi 

A general relaxation time r out was chosen to ac- 
count for the phonon bottleneck effect outside the region 
around H z = 0. The time evolution of the magnetiza- 



tion is given by M(t) = - [M(t) 
for small values of T yields 



M cq (H z (t))] It that 



M(t) = M cq + (M sat - M oq )e 



-t/r 



(11) 



where M sat = M(i? sat ) and M oq = M cq (H z ) denote the 
magnetizations at saturation and thermal equilibrium, 
respectively. H z (t) = H sat + Tt was swept from H z = 
H sat = -10.0 kG to +10.0 kG. r = r in for \H Z \ < 1.5 kG 
and t = T out for \H Z \ > 1.5 kG. The value of H z that 
marks the border between the two regimes is defined by 
the magnitude of H e g tX . 

The butterfly hysteresis in Fig. [3] was best reproduced 
with ri n = 2100 s, whose temperature dependence turned 
out to be insignificant below 10 K. Assuming a realistic 
sound velocity of v = 1400 m/s, 2 » this value corresponds 
to go = 0.72 mK, i.e. about 3 orders of magnitude smaller 
than the values reported for MnizS and Feg££. This low 
value for go in Fe2 arises from the fact that the num- 
ber of thermal phonons is very small due to the weak 
phonon-bath coupling in our experiments. Hence, it does 
not necessarily reflect a weak spin-phonon interaction. 
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Tout = 550 s as obtained from the fit is in overall reason- 
able agreement with the experimental values for T2 de- 
duced from the magnetization relaxation curves (Fig. [3J . 
Importantly, T; n > T out , which leads to the typical fea- 
tures of a butterfly hysteresis loop. The resonant phonon 
field is H s,x = 1 kG, which is about three times larger 
than the value obtained from the analysis of the data in 
Fig. [S] indicating that H c g, x depends on the sweeping 
rate. The best fit was obtained with sweeping rates of 
2.5 G/s for \H X \ < 5 kG and 3.3 G/s for \H Z \ > 5 kG, i.e. 
in very good agreement with the experimental values^ 
To fit the butterfly hystereses at 3 and 5 K (Fig. 0}, 
r out and H c ff tX were multiplied with the Boltzmann fac- 
tor e A(i/fcTi-i/fcT 2 ) j where Ta _ Lg k and Ti = 3 K or 5 
K. The decrease of r ou t and H c s^ x with increasing T re- 
flects the fact that (i) the coherence and thus relevance of 
the resonant phonons decreases, and (ii) the phonon-bath 
coupling increases. 



D. Intrinsic Spin-Phonon Dynamics 

Insight into the intrinsic spin-phonon relaxation prop- 
erties of Fe2 unaffected by the PB is provided by the 
AC susceptibility experiments on samples with a good 
thermal contact with the bath. The absence of a PB sit- 
uation is confirmed by the identical relaxation rates 1/r 
found for single-crystalline (sample B) and powder sam- 
ples (sample C) at given values of Hue, as otherwise 1/r 
would increase with the size of the crystallites as phonon 
scattering at the boundary walls decreases 

With an 5 = 5 ground state and an energy barrier 
of A = 2.40 K for spin reversal, Fe2 features the princi- 
pal requirements necessary for exhibiting SMM behavior. 
However, from the absence of a %" signal at Hue =0we 
conclude that Fe2 does not behave as an SMM at T > 1.8 
K. Obviously at these temperatures the relaxation over 
the small barrier is beyond the detectable range; presum- 
ably the blocking temperature for SMM relaxation lies in 
the mK region. However, by applying a static field Hue, 
strong x" features are readily observed. We ascribe this 
difference to the fact that at Hue 7^ 0, phonons absorbed 
on one side of the barrier can no longer resonate with 
those emitted on the other side, rendering the relaxation 
considerably less efficient. 

The relaxation rate 1/r is temperature dependent, in- 
dicating that the relaxation occurs via a thermal activa- 
tion process. At Hue — 1 kG, the temperature depen- 
dence of 1/r follows the Arrhenius law (Fig. EJ, thus sug- 
gesting an Orbach processiSi This behavior is typically 
observed for SMMs, and often the value for AE is re- 
garded as a lower limit for the SMM barrier height (note 
that x" features at Hue 7^ have been reported for a 
few SMMs including Mn 4f £ V 4 ^ and Mn 12 ^^ for V 4 , 
the x" signal was observed only at i?DC 7^ 0). However, 
in Fe2 the kinetic energy barrier of AE = 12.6 ±0.2 K 
is 4 times larger than the thermodynamic barrier of the 
5 = 5 state at 1 kG, A = 3.1 K (Fig. 0. Importantly, 




FIG. 9: Double well potential for a bias field of i?DC = 1 kG. 
The arrows indicate individual AS = AM ± 1 spin-phonon 
transitions within the combined 5 = 5, 5 = 4 spin system 
and illustrate a likely scenario for spin-phonon relaxation. 



this value lies close to the exchange splitting between the 
5 = 5 and the 5 = 4 states of only 11.2 K 14 (Fig.EJshows 
the anisotropy and Zeeman splittings of these states at 
Hue — 1 kG). This similarity hints toward a relaxation 
process involving transitions to the 5 = 4 state (Fig.[5J). 
While the time-averaged level populations of the spin sys- 
tem are defined by the Boltzmann statistics, the oscillat- 
ing field induces steady-state spin-phonon transitions on 
both sides of the barrier in an attempt to achieve thermal 
equilibrium. At T max ss 4 K, the 5 = 4 state is populated 
by about 5%, and so is the phonon energy spectrum. As- 
suming a Debye model, at 4 K the density of phonon 
states suitable for, e.g., an \S,M) = |5, — 5) — ► |4, — 4) 
transition is about two orders of magnitude higher than 
that for a |5, — 5} — ► |5, — 4} transition, resulting in a 
strongly increased number of available thermal phonons 
for spin-phonon transitions. Also, at 4 K all the 5 = 5 
levels are highly populated. Consequently, a relaxation 
process involving A5 = ±1 rather than A5 = transi- 
tions (Fig.El) becomes favorable. Therefore we postulate 
that the 5 = 5 and 5 = 4 states form a combined sys- 
tem in which the relaxation occurs from the 5 = 5 state 
over the barrier of the 5 = 4 state. The effective energy 
barrier is then 13.8 K and 12.5 K for transitions from 
the deeper to the shallower well and vice versa, respec- 
tively, in good agreement with the kinetic energy barrier 
AE = 12.6±0.2 K. This agreement further indicates that 
A5 = transitions - although still possible - do not play 
a dominant role in the relaxation process. 
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FIG. 10: Argand or Cole-Cole plots of x" vs x' at H DC = 1 
(a) and 5 kG (b) at fixed temperatures (*, 3.0 K; a, 4.0 K; 
o, 5.0 K). For clarity the data are shifted along the x axis. 
For a single relaxation process the data points should lie on 
a semicircle (solid lines) at each temperature. 



Fig. HOfa , shows the so-called Cole-Cole or Argand plot 
for the AC susceptibility data at -£?dc = 1 kG, in which 
X" is plotted vs x' f° r an frequencies at fixed temper- 
atures as indicated. Importantly, for a given temper- 
ature the data do not lie on a semicircle but deviate 
widely^ thus unambiguously indicating the presence of 
a multitude of relaxation processes with different rates. 
These rates can be calculated using Pauli's master equa- 
tion in conjunction with Fermi's golden rule, and are 
obtained as the eigenvalues of the matrix of transition 



rates 



23,34,35,36.37 



The smallest nonzero eigenrate corre- 
sponds to the rate at which the relative populations of 
both wells reach mutual equilibrium, i.e. to the relax- 
ation of the spin system over the barrier. All the faster 
rates are associated with transitions between levels in- 
side each of the two wells. In Mni 2 , e.g., the smallest 
relaxation rate is about four orders of magnitude smaller 
than the second-smallest^ Consequently, Mni2 relaxes 
with a single observable relaxation rate. Obviously, for 
Fe2 the situation is different. The small energetic spread 
of individual spin-phonon transitions leads to a very nar- 
row distribution of eigenrates that are not resolved ex- 
perimentally. Therefore, the relaxation phenomenon ob- 
served in the AC susceptibility measurements of Fe2 does 
not exclusively correspond to the relaxation over the 
5 = 5/5 = 4 effective energy barrier but also reflects 



individual AS ± 1 spin-phonon transitions on either side 
of this barrier. 

At i?DC = 5 kG the situation is distinctly different. 
The 1/t k T dependence (Fig. |JJ is consistent with a 
direct rather than an Orbach process, with a spin-phonon 
transition energy Too < T = 3 — 5 Ki2i Also, the Argand 
plot (Fig. IlUb ) indicates that 1/t essentially consists of 
a single relaxation rate. These results suggest that at 5 
kG, the observed relaxation is governed by the transitions 
between two levels only, which, however, are difficult to 
identify. Nevertheless, a relaxation process involving the 
5 = 4 excited state appears highly unlikely, because (i) 
at 5 kG the 5 = 4 state has lost its barrier entirely, and 
(ii) < T = 3 - 5K. Consequently, the monitored 
relaxation is governed by transitions between two levels 
of the 5 = 5 state. 

The observed increase of x" between 1 and 5 kG 
(Fig. HJ1 Sec. IIII CI) indicates that the relaxation process 
becomes increasingly inefficient with i?DC- A number 
of reasons may contribute to this behavior. Perhaps 
the most important one is that the number of possible 
spin-phonon transitions involved in the relaxation pro- 
cess shrinks markedly with increasing asymmetry of the 
potential well, concomitant with a change from an Or- 
bach to a direct relaxation process. 



V. SUMMARY AND CONCLUSIONS 

This study provides significant insight into the unusual 
dynamics of the magnetic properties of the spin cluster 
[Fe 2 F 9 ] 3 -. 

First, the observed slow magnetization relaxation in 
Fe2 arises from poor thermal contact between the sample 
and the bath. The resulting phonon bottleneck situation 
allows for the indirect observation of quantum tunneling 
of the magnetization through the energy barrier at 1.8 K. 
Importantly, this result illustrates that a PB situation 
can provide insight into the dynamics of the spin sys- 
tem at temperatures far above the blocking temperature 
below which spin-phonon relaxation becomes directly de- 
tectable. 

Second, in the high-spin system Fe2 the hysteresis 
caused by the PB effect adopts a butterfly shape that 
is phcnomenologically analogous to that observed for the 
low-spin systems Vj.5 and Feg. This result indicates that 
despite the energy barrier in Fe2, the spin reversal at 
H z = occurs much faster than the re-equilibration of 
the spin+phonons system with the bath. As QTM is 
rather inefficient in Fe2, it further suggests that ther- 
mal activation is competitive in achieving spin rever- 
sal. Our microscopic model based on the rapid absorp- 
tion/emission of resonant phonons allows for an accurate 
reproduction of the hysteresis curves observed for Fe2- 

Third, AC susceptibility experiments unobstructed by 
the phonon bottleneck allow for direct insight into spin- 
phonon dynamics. Even on the timescale of our AC ex- 
periments, Fe2 does not behave as an SMM at T > 1.8 K. 
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At Hdq = 1 kG an out-of-phase susceptibility is detected 
with a frequency dependence closely resembling the one 
typically observed for SMMs; hovever, the situation in 
Fe2 is substantially more complex and highly unusual. 
The 5 = 5 ground and the 5 = 4 first excited states 
form a combined system in which the relaxation occurs 
over the barrier of the combined system, favored by the 
much higher density of states for thermal phonons suit- 
able for AS = ±1 than AS = transitions. Importantly, 
this result demonstrates that the measured kinetic energy 
barrier can actually be larger than the thermodynamic 
barrier of the ground state if the excited state becomes 
thermally accessible. This result provides an additional 
twist for the interpretation of the kinetic energy barrier of 
single-molecule magnets derived from AC susceptibility 
measurements. 
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We define the quantum state in the generalized rotating 
frame by IV'grot) = U The transformation of the left 
and right hand side of Eq. (|A.2|I then yields 

.M) ( TT ^ d\^ ot ) du- 1 , \ 

= Ift I U ^ 1 7T— IVgrot) , 



dt 



dt 



dt 



Heffl-0) = WeffC/" 1 |Vgrot) • (A.3) 

Combining the left and right hand side of Eq. l|A.3j) leads 
to 



iti 



(V-i^j^ot) + i^ grot) ) = h cS u-' i^ ot ) . 

(A.4) 

Multiplying Eq. i|A.4|) by U from the left results in 



ih 



gJV^rot) 

dt 



UH cS U- 



ihu- 



ot 



IV'grot) 



grot/ 



(A.5) 
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APPENDIX 

Eq. J7J is derived from Eq. © as follows. Generalizing 
the rotating wave approximation^ the Hamiltonian in 
Eq. 10 can be approximated by 



K 



eft 



25D h. 



5,4 







16D Jl4, t 3 

h 3A 9D 







h-*.-A 





h-4,,-5 

25D 



gHBHzS z , 



(Al) 



M+1,M 



M)(S + M + 1) 
To remove the time 



where h 

e iu M+ i,Mt/2 an d liMM+i = h M+hM 
dependence of Hamiltonian (|A.lf) . a unitary transforma- 
tion U is applied to Eq. I|A.1|) . In order to obtain the 
transformation of Ti c s to the generalized rotating frame, 
U has to be applied to the time-dependent Schrodinger 
equation 



at 



(A.2) 



from which we can directly read off the transformed 
Hamiltonian 



dW 



dU, 



Hffr 0t = UHcsU^ - iMJ—- = un cS u^ + ih—U* 1 



dt 



dt 



(A.6) 



We now choose the following unitary transformation: 

5 



u= \ M )( M \ 



(A.7) 



Af=-5 

Evaluating H' s — UHcfiW yields 



25L> h' 5A 
h' 4>5 16D hi <3 
h' 3 4 9D 



'-5,-4 



gHBHzS Z7 



where h\ 



M+l,M 



= h M +i,M^ M ^ M+l)t and h' } 



(A. 



M,M+1 



^m+i,m- Setting u>m+i,m 
the time dependence of 



Wjifli + ujm = eliminates 



n 



grot 



off = K S + ih^tf 



(A.9) 
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